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INTRODUCTION 



The study of atmospheric diffusion can be simply defined as 
the investigation of the diluting effects of the atmosphere on a 
released contaminant. The atmosphere is most often in a 
turbulent state implying that its motions can be predicted only 
in a stochastic sense. Since virtually all gaseous pollutants 
and most airborne particulates very quickly become well mixed on 
the molecular level in the troposphere, diffusion of those 
materials can be described in a similar probabalistic manner. 
Introducing statistics automatically requires the user to 
consider the variations of those statistics in time and space. 

As an example of the later, a predicted crosswind concentration 
distribution represents only a first order approximation to the 
actual distribution, which will often show significant variations 
from that prediction. The study of those fluctuations is 
presently only beginning (see Sawford, 1985). 

This study focuses on the time variations of diffusion 
statistics, specifically, the standard deviation of the crosswind 
concentration distribution (sigma-y). Skupniewicz and Schacher 
(19814b) have shown that overwater releases of material will 
diffuse in the crosswind (y) direction with two scales of motion. 
On one scale, material will disperse about the centerline of the 
plume or cloud due to turbulence of length scales close to the 
plume or cloud size. This scale of dispersion is referred to as 
relative dispersion, implying that the statistics can be observed 
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by the motions of parcels relative to each other in a moving 
frame of reference. On the second and much larger scale, the 
plume or puff will disperse relative to a fixed axis, usually 
chosen to be the mean wind direction. Turbulence acts on the 
instantaneous plume or puff as if it were a single entity, and is 
referred to as single-particle diffusion. A more common and 
understandable synonym is meander. 

Skupniewicz and Schacher (1984a) implicitly considered the 
combined diffusion properties of both components when hourly 
averaged plume parameterizations were derived from tracer 
experiments. In this report, the two components are separated in 
order to supply more information on the nature of the diffusion 
processes, and resultantly more knowledge of potential hazard 
from a contaminant release. 
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2. EXPANDED DATA BASE 



Most descriptions of plume parameters rely on empirical 
formulae and/or major assumptions about the physical behavior of 
the atmosphere. Verification of these parameterizat ions 
therefore should utilize data gathered under a wide variety of 
atmospheric conditions and geographical locations. The 
Environmental Physics Group (EPG) at the Naval Postgraduate 
School ( NPS) has made a significant effort to collect a diverse, 
overwater data base over a several year period for this purpose. 
These data are reported in Skupniewicz and Schacher (1984a) and 
will be hereafter referred to as the NPS data. The NPS data 
primarily consists of concentration profiles collected from an 
aircraft platform at various distances from a continuous surface 
release of SF6 gas. Part of the work reported here was to 
integrate a new data set collected by the German Military 
Geophysical Office (GMGO) (see Groll et al, 1984) into our data 
base. These data were obtained in the North Sea via continuous 
SF6 releases from a ship and subsequent downwind plume transects 
with a second ship sampling gas concentrations. Table 1 lists 
information on the three experiments made available to EPG. 
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Experiment Data No. of Sampling Boats Ho. of Plume Transects 

dAc* A i'i «. i 

23 May-06 Sep 79 3 516 

14 Apr-29 Apr 80 4 558 

i . r. rc y r e- 8*it em;.. • 3c ertoi^qi .jgs./ oe M 

04 Nov-13 Nov 80 4 260 

To ,v o i vene.. x . y i , -j u ? i iub gf 'i J 'i -cri> » :* • *jfr . i ' j 

Table 1. GMGO overwater tracer experiments- 'aVailables^&t- EPG-NPS. 
Each transect produces an "instantaneous" sigma-y value used in 
thi's analysis. ^ . * 



plume 



Table 



Task I specifically involved several steps needed to produce 
parameters from the GMGO data: 

1) transferring original 9-track data files to mass storage 
files in IBM 3033 readable format, 

2) combining half-hour average meteorological files with 
nearly coincident tracer profiles, and 

3) calculating the second moment of the concentration 
distribution (sigma-y) for each profile and the 
downwind distance from navigational information. 

2 lists the resulting data set contents and format. 
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VARIABLE 



DEFINITION 



UNITS 



FORMAT 



DATE 


day, month, year 


no units 


312 


TIME 


hour 


Z 


14 


RH 


relative humidity 


% 


F5. 1 


WD 


wind direction 


deg 


14 


DB1 7 


dry bulb at 17m 


C 


F5. 1 


WSM 


wind speed (MET) 


m/s 


F5.2 


WB1 7 


wet bulb at 17m 


C 


F5. 1 


ET 


water temperature 


C 


F5. 1 


DB3 


dry bulb at 3m 


C 


F5. 1 


WB3 


wet bulb at 3m 


C 


F5.1 


GRAD 


gradient 


C/1 000m 


14 


SB 


speed of boat 


m/s 


F5.2 


SS 


source strength 


m^/hr 


F2.0 


SH 


sampling height 


m 


12 


w ss 


wind speed (SF6) 


m/s 


F5.2 


DS 


distance from source 


m 


F8. 1 


BTS 


bearing toward source 


deg 


F5.0 


HD 


head i ng 


deg 


F2.0 


# pt 


no. of points used for CWCI 


no units 


12 


CWCI 


cross wind concentration 
integration 


ppb-m 


F8.0 


MEAN 


mean mass position 


m 


F6. 1 


SD 


standard deviation 


m 


F6 . 1 


Table 2. 


Contents of new GMGO overwater 


tracer data set 


added to 


the EPG overwater data base. Format is 
mass storage data set name is "MSS.F3896 


in FORTRAN code. 
. DATA. NEWSM". 


Public 
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3. RELATIVE DIFFUSION PARAMETERIZATION 



The basic approach used in parameterizing the relative 
diffusion is to group observations according to Pasquill-Gif ford 
equivalent stability classes (see Turner, 1967), and then regress 
sigma-y versus range (from the release point) for each class. 
Sigma-y was defined in the last section as the second moment of 
the concentration distribution for a given profile. The 
customary assumption that concentration in the y direction is 
independent of concentration in either the x or z direction leads 
to the conclusion that these instantaneous "snapshots" of the SF6 
plume are identical to the crosswind dimensions of a hypothetical 
"puff" or "burst" release under the same atmospheric conditions. 

The methodology for determining stability class is described 
in detail in Skupniewicz and Schacher ( 1 98-Ma) . Briefly, the 
scheme uses windspeed, air-sea temperature differences and 
relative humidity to produce stability classes equivalent to the 
well-known Pasquill-Gif ford diffusion categories. [Whether a true 
equivalence exists is a matter of controversy and is discussed in 
Skupniewicz and Schacher (198-4a) . 3 

Under this scheme classes are assigned a letter designating 
atmospheric stability, with "B" representing the most unstable 
situations, "C" representing moderately unstable, "D" representing 
neutral, and "E" representing moderately stable conditions. 

After dividing the data according to stability, the subsets 
were examined for dependence upon wind speed, time of day, and 
day to day trends. No obvious dependences were found. Figures 
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1-3 show sigma-y vs. downwind distance, with the integer 
truncation of the wind speed used as the point designator. 

In as much as windspeed gives a first order approximation to 
sea-state, increased surface stress within a given stability 
class has a negligible effect on diffusion. This result is 
particularly surprising for the neutral case (class D) , where a 
small air-sea temperature difference will always produce neutral 
stability under this scheme regardless of windspeed. Intuitively, 
one would expect that higher windspeed would enhance the plume 
spread. Because stable or unstable profiles tend towards 
neutrality as windspeed increases in this stability scheme, 
windspeed dependence is not expected in non-neutral categories. 

It is also noteworthy that stable situations are found 
exclusively in light wind situations in the North Sea experiments. 
These situations were not found in the California experiments. 

Based on the above data, the parameterization chosen for 
relative diffusion is a simple linear relationship. The 
following factors influenced this choice. Instantaneous releases 
into a turbulent atmosphere will theoretically experience the 
following four different growth regimes: 



a ) 


the 


initial stage 


b) 


the 


inertial stage 


c ) 


the 


central stage 


d) 


the 


final stage 
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Figure 1, Sigma-y relative vs, downwind distance. Symbol of 
data point is integer truncation of windspeed in m/s. Only 
data 3200 m from source is plotted. 

a) Pasquill-Gif f ord stability class B (moderately unstable), 
GMGO data. 23 data points are hidden. 

b) Pasquill-Gif ford stability class C (slightly unstable) , 
GMGO data. 15 data points hidden, 1 data point missing. 
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Figure 2. Same as figure 1, except a) Pasquill-Gif f ord stability 
class D (neutral), GMGO data. 194 data points hidden, 114 data 
points missing. b) class E (slightly stable), GMGO data. 8 data 
points hidden, 153 data points missing. 
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Figure 3. Same as figure 1, except a) Pasquill-Gif f ord stability 
class D (neutral), NPS data. 13 data points hidden, 145 data points 
missing. b) class E (slightly stable), NPS data. 1 data point 
hidden, 45 data points missing. 
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The initial stage applies to distances of only 10-IOOm. The 
inertial stage will be relevent from 100m to less than 1km. The 
central stage will apply from 1km to 10km and the final stage 
will apply thereafter. These bounds are highly variable. 

The initial stage of growth is valid when the instantaneous 
puff "remembers" its initial size. In this regime cloud size 
will increase linearly with distance or travel time because each 
particle in the cloud will move with its initial velocity and the 
particles have not had time enough to become correlated with each 
other . 

The initial stage is followed by growth dominated by 
turbulence in the inertial subrange. In the inertial region, 
plume growth can be shown to follow 

0y r = aet 3 , (1 ) 

where 0 y r is the relative sigma-y, 
a is a constant, 

e is the turbulent kinetic energy dissipation rate, and 
t is travel time of the cloud. 

This inertial stage of growth is valid only where the width of 
the cloud is small compared to the average height of the cloud; 
surface-released clouds seldom meet this criterium. 
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Smith and Hay (1961) laid the foundation for considering 
cloud spread as a function of turbulent energy in a "sliding" 
spectral window. They found both experimentally and 
theoretically that cloud spread will approximately follow a 
linear growth relationship 

<r =0 . 22 ix , ( 2 ) 

yr 



where i is the total turbulence intensity, and 

• = 1/2 
1 0 " 



(3) 



where x is downwind distance, 

u,v,w are fluctuating wind components, and 
U is mean windspeed. 



Mikkelsen and Eckman (1 
results from overland e 
to medium ranges. 

At long distances 
hypothetical particles 
spread of the cloud wil 
particle solution, 



98^1) have recently given support to these 
xperiments for surface releases in short 

from the source, correlation between 
within the cloud approaches zero, and the 
1 behave like the asymptotic single- 



0 « t 1/2 

y 



( 4 ) 
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Large scale motions, however, are often organized into coherent 
vortex structures (i.e. cyclones). As a result, this limit is 
rarely reached, and growth continues along a more or less linear 
asymptote . 

All of these above arguments and the rather linear shape of 
the data scatter lead us to a linear paramterization. The 
results are shown in figures JJ-6 and summarized in table 3. Note 
that there is very little unstable data at ranges greater than 
3000m. This is primarily due to the fact that, under these 
conditions, diffusion is enhanced and tracer concentrations are 
greatly reduced far from the source. This makes the plume 
difficult to locate, and may bias measurements towards narrow, 
concentrated profiles. The most striking feature is that there 
is no clear distinction between P-G classes B and C ( GMGO data), 
or between classes D and E (both data sets). Those pairs are so 
closely matched that we suggest only a two-class parameterization 
(also shown in table 3) ; one for all unstable conditions and a 
second for neutral to stable stratification. While turbulence 
intensity is not explicitly used in the paramterization, this 
two-class system is supported by the calculated wind variances 
from the NPS data shown in table 4. These wind data include some 
measurements from periods not coincident with tracer releases, 
and do not include some periods when tracers were released. 
Nonetheless, the largest change in turbulence (when moving from 
one to an adjacent category) occurs between class D and class C. 
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Figure 4. Sigma-y relative vs. downwind distance with solid line 
representing the recommended parameterization of table 3. Symbol 
A represents one data point, B represents two data points, etc. 

Same data as figure 1, except long range data is included. Data 
is a) Pasquill-Gif ford stability class B, GMGO data, and b) class 
C, GMGO data. 
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Figure 5. Same as figure 4, 
class D, GMGO data, and b) 



except a) Pasquill-Gif f ord stability 
class E, GMGO data. 
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Figure 6. Same as figure 4, except 
class D, NPS data, and b) class E, 



a) Pasquill-Gif ford stability 
NPS data. 
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DATA 


PASQUILL-GIFFORD CLASS SLOPE 


//POINTS 


GMGO 




B .01167 


1 32 


GMGO 




C .0357 


1 110 


GMGO 




D .0227 


603 


GMGO 




E .0178 


206 


NPS 




D .0178 


215 


NPS 




E .0228 


72 


GMGO 




B + C ** .01110 


272 


GMGO + NP 


s 


D+E ** .0211 


1 096 


**Last 2 


rows are 


suggested values. 




Table 3. 


Relat i ve 


diffusion parameterization for 


overwater , 


surface , 


medium range (1-12km) releases. Slopes 


are least 


squared 1 


inear regression to data. Sigma-y relative = 


( Slope )x ( 


downwind 


distance ) . 





Pasquill-Gifford 
Stability Class 


Sigma Wind 
1 min ave 


Direction (deg) 
10 min ave 


// Hours 
Of Data 


B 


3-96 


6.41 


7 


C 


2.95 


H.59 


1 1 


D 


1 .82 


2.92 


1 1 it 


E 


1 .32 


2.61 


1 1 



Table 4. Sigma theta 
Sampling rate is 1 Hz 



as a function of P-G stability class. 
Data is from NPS experiments. 



This two-class feature can also qualitively seen in figure 7 
where turbulence intensity as a function of stability (defined 
using the Mon i n-Obuknov length [ L ] ) was plotted by Schacher et 
al (1982). These data show a step increase in sigma theta 
(turbulence intensity) for the 1 minute averages at roughly L = 
-40m. 1 hour averages are also plotted for a comparison, even 

though 1 minute averages are much more applicable to relative 
diffusion in these experiments. Note that the 1 hour values do 
not show this effect and are 4 to 5 times larger than 1 minute 
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values . 



This suggests that our two-class approach may break down 
as the cloud travel time increases. We therefore do not recommend 
interpolating our results to distances greater than 10 km. 
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Figure 7. Wind direction sigma (sigma theta) vs. z/L 
(z is the height, 10m, and L is the Monin- 
Obukhov length) . X's are 1 min averages, 
and O's are 1 hr averages. Data is from 
HPS. (see Schacher, et al. , 1932) 
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4 . MEANDER (SINGLE PARTICLE) PARAMETERIZATION 



For purposes of describing meander, information on the motion 
of the plume is available only from the NPS data set. Meander was 
measured by calculating the distance between the center of mass 
for a given profile and the axis defined by the mean wind vector. 
The method used to calculate the mean wind vector was crucial. 
After experimenting with various techniques, we decided to let 
the average travel time of the plume define the proper averaging 
period for the vector; this travel time was roughly one-half hour. 
Since the wind was measured only at the release point, tracer 
profiles were correlated with average wind vectors in half-hour 
bins offset by 15 minutes to account for plume transport. 

The NPS experiment was designed so that the plume was sampled 
at nearly discrete ranges from the release point. This sampling 
procedure allowed us to segregate data into range bins. The 
collection of off axis center of mass positions for a given range 
bin represents a probability distribution for the meander of a 
plume or puff. Each experimental day produced several distribu- 
tions at various downwind distances from the source. The distri- 
butions were then analyzed for normality, and the standard devia- 
tions (which represent sigma-y due to meander) were determined. 

The data were examined for dependence on windspeed, 
stability class, inversion height, experimental day, and time of 
day. The day to day variability was so large that it tended to 
swamp other dependencies. Figures 8-21 show day by day plots of 
the center of mass distributions interlaced with the respective 
wind direction time histories and spectral plots. 
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DOWNWIND DISTANCE (m) 



DEVIATION FROM AVERAGE HIND DIRECTION (DEG) LOG (nS (n) ) 





Figure 9. a) Wind direction tine series, and b) pcwer spectra 
for tracer release depicted in figure 8. Time series data are 
2 minute averages of wind direction sairpled at 1 Hz. Spectral 
points are smoothed with a 5% sliding average. S(n) is wind dir- 
ection power spectral density at frequency n. D is the total 
sanpling time. < DIR? is the mean wind direction. 
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Figure 10. Same as figure 8, except 22 June 1982. 
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Figure 11. a) Wind direction time series, andb) power spectra 
for tracer release depicted in figure 10. 
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io0 | Figure 12. Same as figure 8, except 24 June 1982. 

| 10 observations hidden. 
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Figure 13. a) Wind direction time series, and b) power spectra 
for tracer release depicted in figure 12. 
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Figure 15. a) Wind direction time series, and b) power spectra 
for tracer release depicted in figure 14. 
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Figure 16. Sane as figure 8, except 27 June 1982. 2 observations hiddsn. 
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Figure 17. a) Wind direction time series, and b) power spectra 
for tracer release depicted in figure 16. 
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Figure 18. Sarte as figure 8, except 28 June 1982. 
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Figure 19. a) Wind direction time series, and b) power spectra 
for tracer release depicted in figure 18. 
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Figure 20. Sane as figure 8, except 29 June 1982. 
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Figure 21. a) Wind direction time series, and b) power spectra 
for tracer release depicted in figure 20. 
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Individual data points are marked with integer values of the 
measurement hour for a given experimental day. Note that the 
mean position of the cluster often moves off the centerline, but 
data are distributed evenly over the course of a day. Large 
variations in the meander "envelope” can be seen from day to day, 
and the wind direction time series differ radically. The power 
spectra also vary greatly from day to day. Spectral parameters 
are summarized in table 5. Spectral gaps (defined in the table 5 
caption) exist in most situations, with slopes of approximately 
-2 in the low frequency part of the spectra. 

There appears to be no clear relationship between the total 
energy in the turbulent and low frequency parts of the spectra. 
Since meander is driven by the low frequency energy, this 
discourages attempts at parameter izations of meander based on 
surface layer scaling. 

As stated above, center of mass distributions were grouped 
by experimental day and range bin. Variations of the position of 
the measurement aircraft dictated range bins to be 2 km long. 
Relevant statistics are given in table 6. Examination of the 
Shapiro-Wilk statistics (Shapiro and Wilk, 1965) indicates that 
the crosswind distributions are represented fairly well by normal 
distributions. This suggests that Gaussian models will be valid 
even for dispersion caused by meander, given the proper sigma-y. 
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Dat 


e 


/21 


/82 


/22 


/82 


/2H 


/82 


/25 


/82 


/27 


oo 

r\> 


/28 


/82 


/29 


/82 



Classicali peak 
10 min 
20 
30 
20 
5 

10 



Ga £2. 


Slopes 


20 min 


-6.5/3 


50 


-7/3 


50 


-6.5/3 


25 


-6/3 


15 


-5.5/3 


15 


-5.5/3 



? ? 



1 Classical peak is the inverse of the frequency associated with 
the "turbulent energy producing" subrange maximum at the low 
frequency end of the -5/3 slope (inertial subrange) region. 

2 Gap is the inverse of the frequency where a minimum occurs on 
the low frequency side of the "turbulent energy producing" 
subrange . 

3 Slope is the log (power) vs. log (frequency) relationship in 
only the "low frequency" part of the spectrum. Gap 
frequencies are not included. 



Table 5. Summary of wind 
in figures 9b, 11b, 13 b, 



direction spectra 
5b , 17b, 19b, and 



for NPS data depicted 
21 b. 
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DATE RANGE(km) N 


MEAN! (jn) 


STD DEV? ( 


ml SKEW? 


KURTOSISn W^ 


21 Jun82 


2 


25 


-293.0 


1*19.7 


-0.2861 


-0.3701 


0.9705 


21 Jun82 


*1 


26 


-356.5 


255.6 


+0.6759 


+0.0095 


0.955*1 


22 Jun82 


*1 


*13 


-1733-7 


392.1 


-0.0101 


-0.68*15 


0.9721 


2*1 Jun82 


2 


13 


-80.5 


110.3 


+0. 88*10 


+0.6898 


0.9077 


2*JJun82 


*1 


3*1 


-678.0 


391 .5 


-0. *1567 


+ 1 .0885 


0.9590 


25 Jun82 


2 


1 *J 


-201 .3 


50. *1 


+0.733 1 * 


-0.82*12 


0.8912 


25 Jun82 


*1 


38 


-7*10. 9 


166.3 


+0.7936 


+0.8523 


0 . 9 * 17*1 


27 Jun82 


2 


9 


-185.9 


103.3 


-0.0621 


-0.0*1*11 


0.9567 


27Jun82 


*1 


17 


-895. *1 


198.8 


-0. 1 5*11 


-0.122*1 


0.967*1 


27 Jun82 


6 


*1 


-1516.6 


*11 5. *1 


+ 1 .22*12 


-0.7663 


0.8893 


28Jun82 


2 


6 


-*117.3 


83. *1 


+0.8019 


- 0.1 838 


0.9396 


28 Jun82 


*1 


12 


-795.0 


260.9 


-0.1 653 


+ 1 .*1008 


0 . 9*170 


28 Jun82 


6 


6 


-970.6 


160.5 


-1 . 3268 


+ 1 .7325 


0.8868 


28 Jun82 


8 


5 


-1165.7 


*167.0 


+ 2.1 2*18 


+ *J . 60*13 


0.6790 


29 Jun82 


2 


3 


-613.9 


*177.9 


+0.9695 




0.9613 


29 Jun82 


*1 


12 


-913.0 


1008.9 


+0.6763 


+0.0379 


0.937*1 


29 Jun82 


6 


9 


-789.9 


966.7 


+0. 8*135 


-1 .0370 


0.8362 


29 Jun82 


8 


1 1 


-1177.6 


1380.2 


+ 1 .21 22 


+0.7888 


0.8509 



1 MEAN: Mean position of plume center (as measured from the 

mean wind centerline) of mass distribution for a 
given range bin. 



2 STD DEV: Standard deviation of plume center of mass for a 

given range bin. 

3 SKEWNESS: Measures the tendency of deviations to be larger in 

one direction than the other. Values are unbounded. 
Positive values indicate a "tail" of values exists 
in the values larger than the mean, negative values 
indicate a "tail" in smaller values. 



*1 KURTOSIS: Measures the heaviness of the distribution "tails" 

(-2.0 < Kurtosis < °°). Relatively large values 
indicate normality assumption should be questioned, 
but this statistic is unstable for a small number 
of samples which is characteristic of this 
experiment . 



5 W: Shapiro-Wilk statistic for the center of mass 

distribution. Measures the normality of the 
distribution (0.0 < W < 1.0). Values close to 1.0 
indicate more normal distributions. 

Table 6. Meander statistics for NPS meander data . N is the 
number of profiles for a given crosswind center of mass 
d i str ibut ion . 
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While we were not able to parameterize meander with mean 



surface layer meteorological variables, direct measurements of 
lateral turbulence intensity did correlate well with the observed 
meander standard deviations. Taylor (1921) predicts that when 
the length scale of turbulence is much larger than the scale of 
the plume (sigma-y meander in this case) the following simple 
formula will apply: 



o 

ym 



i x 
v 



(5) 



where 0 y m is the single-particle lateral standard deviation, or 
sigma-y meander, 

i v is the lateral turbulence intensity, 
x is downwind distance. 

This is often referred to as the "near-field" approximation to 
single-particle diffusion. For application to the present 
analysis, we define the lateral turbulence intensity as 



„ 2 S ‘4 

<v > 2 



( 6 ) 



where <v 2 > is the variance of the crosswind component measured 

over the travel time of the tracer (approximately 1/2 
hour) and averaged over the entire day. 

U is the mean wind speed 
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Figures 22-28 show equation (5) plotted for each 
experimental day with the standard deviations of table 5 (sigma-y 
meander) depicted as data points. Also shown is a least-squared 
linear fit to the data. 

Note the good agreement between this simple model and the 
observed data in almost every case. The key to success is 
matching the travel time of the plume to the RMS window in the 
turbulence intensity measurement. Also, since the measured 
variance is highly non-stationary , it is necessary to average 
this quantity over a full day. 

The agreement between the data and the near-field 
approximation is also supported by our spectral results. 

The spectra continue to increase in power through the entire 
spectral window (4-8 hr ) which suggests that turbulent motions 
operating at these very low frequencies have very long time 
scales as equation (5) assumes. 
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meters °w~ me * ers 




Figure 22. Sigma-y meander vs. downwind distance. Error bars 
are the standard deviations of downwind distance positions for 
a given range bin. Upper number is Shapiro-Wilk statistic and 
lower is number of data points. Solid line is best fit to data. 
Dashed line is equation 5. Data is NPS, 21 June 1982. 




Figure 23. Same as figure 22, except NPS data, 22 June 1982. 
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* meters <7 - meters V melers 




DWD - meters 

Figure 24. Same as figure 22, except NPS data, 24 June 1982. 




DWD - meters 

Figure 25. Same as figure 22, except NPS data, 25 June 1982. 
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Figure 27. Same as figure 22, except NPS data, 28 June 1982. 




Figure 28. Same as figure 22, except NPS data, 29 June 1982. 
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5. CHEMICAL WEAPONS HAZARD FORECAST PROGRAM MODIFICATIONS 



Part of the current research work is to modify the Chemical 
Weapons Hazard Forecast Program (CWHFP). This section assumes 
the reader is familiar with that program and understands the 
program's design, inputs, and operation. EPG modified the 24 
May 1984 "test and evaluation" version of this program supplied 
to EPG by J. Branum (1984). Later versions of this model may 
exist, and the following modifications can be easily transplanted 
into such hybrids. 

The proposed modifications specifically make use of the 
parameterizations described above and are designed so as to 
supply the user with a better understanding of the character of 
atmospheric diffusion processes, and hopefully, more information 
on how to avoid hazards. The previous model used meteorological 
inputs to calculate a stability category, and then calculate one 
hour averaged plume dimensions based on stability dependent 
sigma-y and sigma-z algorithms. This version uses the same 
stability scheme, but replaces the one-hour sigma-y functions 
with the relative diffusion parameterization presented above. 
Sigma-z parameterizations are unchanged. The resultant 
isopleths represent either "dosage" values from an instantaneous 
"burst" release or a concentration "snapshot" from a continuous 
plume . 
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In order to present the meander of a plume (offaxis 
excursions of a puff) the single-particle parameterization 
developed above is used to obtain a "meander envelope." This 
envelope is simply the 2 sigma-y position of the plume or puff 
center of mass distribution represent ing- a 95? probability of 
impact. This envelope is superimposed on the instantaneous 
isopleths in order to show the combined impact of the two 
diffusion processes. An artificial "ripple" is convolved into 
the instantaneous solution to offer the user a visualization of 
meander. The frequency and amplitude of this ripple has no 
physical basis, and is only installed to help the user understand 
the nature of the meander process. 

Figures 29-36 show several examples of model output. Labels 
display lethal dose levels on the instantaneous isopleths and a 
footnote explains the meander envelope. Instruction manuals 
should explain the arbitrary nature of the ripple in order to 
avoid false conclusions based on these plots. Note the change in 
the frequency of the ripple as the meander envelope shrinks (see 
figures 29, 30, 31 or 32, 33, 3*0. This change again was 
arbitrary , and was only intended to suggest the possibility that 
the important eddy sizes decrease as the meander decreases. The 
amplitude of the ripple, again, was arbitrary (selected to be 1 
sigma-y meander) and suggests that the instantaneous position of 
the cloud is somewhat random and rarely reaches the meander 
envelope bounds. Also note the change in the instantaneous cloud 
dimensions with a change in stability class (figures 29, 30, 31, 
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CHEMICAL WEAPON HAZARD FORECAST PROGRAM - SNAP O.a 




TERRAIN TYPE 
MEAN WIND 

SIGMA WIND DIRECTION 
MEAN AIR TEMP 

mean sea-surface TEMP 

RELATIVE HUMIDITY 
SURFACE MIXING LAYER HT 
STABILITY CATAGORY 
MUNITION TYPE 
SOURCE TYPE 

SOURCE SIZE effect i ve) 
SOURCE RATE 



OPEN-SEA 

10 KTS FROM 275 DEG TRUE 

20 DEG 

10 DEG C 

17 DEC C 

70 * 

0 METERS BASED ON •* DEFAULT ** 

B MODIFIED PASOUILL 

M> 116-SIZE BOMB/MISSILE 

POINT-BURST 

88 KG 

INSTANTANEOUS 



BSS3SS5SSSsrSS=C=CCCSS83=CS3SSS=SSS3=CCS:5S33CSCSCS 



3SSSSSSSSS3SSSSSSSCCSSCS 



CONTOUR LABEL POTENTIAL CASUALTY EFFECTS APPROX MAX 

* DOSE-AGENT » ■ WITHOUT PROTECT I ON' RANGE 



- LD50-GP 


50*. 


DEATHS - MOST 


INCAPACITATED 


875 


*j ARDS 


• LD1-GD 


r. 


DEATHS - many 


I NCAPAC I TATED 


1933 


YARDS 


- LDO-GP 


NO 


deaths - SOME 


I NCAPAC ITATED 


2588 


YARDS 



OUTS IDE - THICK* CONTOUR IS 95*. PROBABILITY THAT CENTER OF MASS WILL CROSS 



Figure 29 • Chemical Weapons Hazard Forecast Program output 
sample — unstable atmosphere, large meander* 
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9NAP XX. X 



CHEMICAL WEAPON HAZAPD TOftECRST PROGRAM - 




TERPRIN TYPE 
HERN WIND 

S1GMR WIND DIRECTION 
HERN RIR TEMP 
MERN SER - SUPFRCE TEMP 
RELRTIVE HUMIDITY 
SUPFRCE MIXING LRYEP HT 
STABILITY CRT RGORY 
MUNITION TYPE 
SOURCE TYPE 

SOURCE SIZE Ktf f«ct 1 > 
SOURCE RRTE 



OPEN-SEA 




10 


KTS FROM 


275 DEG TRUE 


ie 


DEG 




to 


DEG C 




1 7 


DEG C 




70 


V. 




0 


METEPS BASED ON *♦ 


defrult t* 


B 


MODIFIED PAS&UILL 




MK 1 


1 E-SIZE BOMB ^MISSILE 




POI 


NT-BURST 




88 


KG 




INS 


TRNTRNEOUS 





SSB3SXSXS3X3S 






CONTOUR' LRBEL POTENTIAL CRSURLTY EFFECTS RPPROX MRX 

<DOSE-RGENT ) 'WITHOUT PROTECTION* RRNGE 



- LD50-GP 


50*. 


DERTHS 


- MOST 


INCRPRC ITRTED 


875 


YRPDS 


- LD1-GD 


1 *. 


DERTHS 


- many 


INCRFRC I TRTED 


1933 


YARDS 


- LDO-GD 


NO 


DERTHS 


- SOME 


INCRPRC ITRTED 


2588 


*i RPDS 



OUTSIDE -THICK • CONTOUR IS 95*. PROBABILITY T HR T CENTER OF MASS WILL CROSS 



Figure 30. Chemical Weapons Hazard Forecast Program output 
sample — unstable atmosphere, moderate meander. 
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CHEMICAL WEAPON HAZARD T0PECR9T PROGRAM - GNfiP XV. X 




TERRAIN TYPE 
MEAN WIND 

SIGMA WIND DIRECTION 
HERN RIP TEMP 
HERN SEA-SURF RCE TEMP 
PELRTIVE HUMIDITY 
SURFACE MIXING LRYEP HT 
STABILITY crtrgory 
MUNITION TYPE 
SOURCE TYPE 

SOURCE SIZE • effective) 
SOURCE RATE 



10 

7 

10 

1 ? 

70 

0 



ITS 

DEG 

DEG 

DEG 



FROM 275 DEG TRUE 



METERS BASED ON ** DEFAULT ** 
B MODIFIED PASOU ILL 
MY 1 16-SIZE BO MB' MISSILE 
POINT-BURST 
8P KG 

instantaneous 






Ksxssesasssscssssa 



a x s = » » * 



asscsxaa 



CONTOUR LABEL 
< DOSE-AGENT > 



POTENTIAL CASUALTY EFFECTS 
i WITHOUT PROTECTION) 



APPROX MAX 
RANGE 



- LD50-GD 

- LD1-GD 

- LD0-GD 



50*. DEATHS - MOST INCAPACITATED 
1*. DEATHS - MANY INCAPACITATED 
NO DEATHS - SOME INCAPACITATED 



875 YARDS 
1933 YARDS 
2588 YARDS 



OUTSIDE 'THICK' CONTOUR IS ?5\ PROBABILITY THAT CENTER OF MASS WILL CROSS 



Figure 31, Chemical Weapons Hazard Forecast Program output 
sample — unstable atmosphere, small meander. 
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CHEMICAL NCR PON HAZARD FORECAST PROGRAM - SNAP XX . Y 




TERRAIN TYPE 
MEAN wind 

SIGMR WIND DIRECTION 
MERN RIR TEMP 
MERN SEA-SURFACE TEMP 
RELATIVE HUMIDITY 
SURFACE MIXING L RYE P HT 
STABILITY CRTRGORY 
MUNITION TYPE 
SOURCE TYPE 

SOURCE SIZE 'effect i > 
SOURCE RRTE 



OPEN-SEA 

10 KTS FROM 27? DEG TRUE 

20 DEC 

17 DEG C 

10 DEG C 

70 X 

0 METERS BRSED ON ** DEFRULT ** 

E MODIFIED PRSOUILL 

MK1U.-SIZE EOMB'MISSILE 

POINT-BURST 

88 KG 

INSTRNTRNEOUS 



XS==S=BSS=SSSSSa=SSS=SS==BZSBS«BSK 






CONTOUR LRBEL POTENT IRL CRSURLTV EFFECTS 

« DOSE-AGENT) 'WITHOUT PROTECTION) 



RPPROX MRX 
RRNCE 



- LD?0-GD 

- LD1-GD 

- LD0-GD 



50’. BERTHS - MOST INCRPRC ITRTED 
1\ BERTHS - MANY INCRPRC I TRTED 
NO BERTHS - SOME INCRPRC I TRTED 



5778 YARDS 
13597 YRRDS 
18638 YARDS 



OUTSIDE 'THICK > CONTOUR IS 95*. PROBABILITY THAT CENTER OF MASS WILL CROSS 



Figure 32. Chemical Weapons Hazard Forecast Program output 
sample — stable atmosphere, large meander. 
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GNPP XX.y 



CHEMICfiL MERPON HRZRPD TOPECRST PROGRAM - 




TERRRIN TYPE 
HERN WIND 

SIGMA WIND DIRECTION 
MERN RIP TEMP 
MERN SER-SURFRCE TEMP 
PELRTIVE HUMIDITY 
SUPFRCE MIXING LRYER HT 

strbility crtrgory 

MUNITION TYPE 
SOURCE TYPE 

SOURCE SIZE l€ff*ctiv»> 
SOURCE PhTE 



OPEN-SER 

10 KTS FROM 273 DEG TRUE 

1 0 DEG 

17 DEC C 

10 DEG C 

70 

0 METERS BRSED ON ** BEFRUIT ♦* 

E MODIFIED PRSOUILL 

MK 1 16-SIZE BOMB-MISSILE 

POINT-BURST 

88 KG 

INSTANTANEOUS 



= = x*»rssc = a** = = = = = = *s« = 5s = = = = = = = = = = = = *s£r»»=:s = »as*sa = a = = »** = = = * = *** = ***a**e**r* 



CONTOUF LR BEL POTENTIRL CASUALTY EFFECTS 

< DOSE-AGENT) .WITHOUT PROTECTION* 



APPROX MR X 
PRNGE 



- LD30-GD 

- LD1-GD 

- LDO-CD 



30*. BERTHS - MOST I NC RPRC I TRTED 
I\ BERTHS - MRNY INCRPRC I TRTED 
HO BERTHS - SOME I NCftPfiC I TRTED 



5776 YARDS 
13597 YAPDS 
18838 YARDS 



OUTSIDE 'THICF' CONTOUR IS 93% FROBRBILIT*. TmRT CENTER OF MRSS WILL CROSS 



Figure 33. Chemical Weapons 
sample — stable atmosphere, 



Hazard Forecast Program output 
moderate meander. 
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CHEMICAL NERPON HRZPPD TOPECflST PPCGPRM - SNfil *a.a 




TERRAIN TYPE 
HERN WIND 

SIGMA WIND DIRECTION 
MEAN AIR TEMP 
MEAN SEA-SURFACE TEMP 
RELATIVE HUMIDITY 
SURFACE MIXING LAYER HT 
STABILITY CATAGORY 
MUNITION TYPE 
SOURCE TYPE 

SOURCE SIZE (effect tut ) 
SOURCE RATE 



OPEN-SEA 

10 KTS FROM 275 DEG TRUE 

7 DEG 

17 DEG C 

10 DEG C 

70 5i 

0 METERS BASED ON ** DEFAULT *♦ 

E MODIFIED PASOUILL 

MK l 16-SIZE BOMB MISSILE 

POINT-BURST 

88 KG 

INSTANTANEOUS 



CONTOUR LABEL POTENTIAL CASUALTY EFFECTS 

‘ BOSE-hGE NTT 'WITHOUT PROTECTION > 



APPROX MAX 
RANGE 



- LD50-GD 

- LB 1 -CD 

- LDO-GD 



50\ DEATHS - MOST INCAPACITATED 
r. DEATHS - MANY I NC APAC I T ATED 
NO DEATHS - SOME INCAPACITATED 



5778 SARDS 
13597 "t ARDS 
18638 YARDS 



OUTSIDE 'THICK) CONTOUR 



95*. 



PROBABILITY that center of 



MASS WILL CROSS 



Figure 34. Chemical Weapons Hazard Forecast Program output 
sample — stable atmosphere, small meander. 
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CHEMICAL UEAPON HAZAPD FORECAST PROGRAM - SNAP XX. V 




TERRAIN TYPE 


OPEN- 


SEA 




MEAN WIND 


5 


KTS 


FROM 


SICMA WIND DIRECTION 


12 


DEG 




MEAN AIR TEMP 


10 


DEG C 




MEAN SEA-SURFACE TEMP 


17 


DEG C 




RELATIVE HUMIDITY 


70 






SURFACE MIXING LAYER HT 


0 


METERS BASED 


ON #* 


STABILITY CATAGORY 


B 


MODIFIED PASQUILL 


MUNITION TYPE 


MK1 16 


-SIZE BOMB'MISSILE 


SOURCE TYPE 


POINT 


-BURST 




SOURCE SIZE '.effective) 


88 


KG 




SOURCE RATE 


INSTANTANEOUS 





275 DEC TRUE 



ssmmssssssmsssssssssssssssmssmssssmmsms 



CONTOUP LABEL POTENTIRL CASUALTY EFFECTS 

(. DOS E-RCEMT ) ■ HI THOUT PROTECT I ON) 



APPROX MAX 
RRNCE 



- LD50-GD 

- L D 1 -CD 

- LD0-CD 



50*. DEPTHS - MOST I NCRPRC I TRTED 
1*. DEPTHS - MR NY INCRPRC I TRTED 
NO DERTHS - SOME I NCRPRC I TRTED 



1301 YRRDS 
2872 YARDS 
3846 YARDS 



OUTSIDE 'THICK* CONTOUR IS 95*. PROBABILITY THAT CENTER OF MASS WILL CROSS 



Figure 35 
sample — 



Chemical 
light wind 



Weapons Hazard Forecast Program output 
conditions . 



51 



CHEMICAL LCRPON HAZftPD TORCCR5T PROGRAM - SNRP XX. X 




TERRAIN TYPE 
HERN WIND 

SIGMA WIND DIRECTION 
MEAN RIP TEMP 
NERN SEA-SURFACE TEMP 
RELATIVE HUMIDITY 
SURFACE MIXING LAYER HT 
STABILITY CRT RGORY 
MUNITION TYPE 
SOURCE TYPE 

SOURCE SIZE <tff»cnvt) 
SOURCE RATE 



OPEN-SER 

20 KTS FROM 275 DEG TRUE 

4 DEG 

17 DEG C 

10 DEG C 

70 * 

0 METERS BASED ON ** DEFAULT ** 

D MODIFIED PR SOU ILL 

MK116-S1ZE BGMB'MISSILE 

POINT-BURST 

86 KG 

INSTANTANEOUS 



3SSBCSSSSSaB*SSSSSSB3*XXSSSSSZSSSSSSS3S2SXS8S3SSSB3S3SSlS133S*SXSttI3SS3S2SXXSSS 



CONTOUR LABEL POTENTIAL CASUALTY EFFECTS APPROX MAX 

' DOSE-AGENT > ‘WITHOUT PROTECTION' PRNGE 



- LD50-GD 


50‘. 


DEATHS - MOST 


INCAPACITATED 


2443 


YARDS 


- lDI-GD 


r. 


DEATHS - MANY 


INCAPACITATED 


5660 


YARDS 


- LD8-CD 


NO 


DEATHS - SOME 


INCAPACITATED 


7714 


YARDS 



OUTSIDE ^THICK' CONTOUR IS 95\ RPOBAE I L 1 TV THAT CENTER OF MASS WILL CROSS 



Figure 36. Chemical Weapons Hazard Forecast Program output 
sample — strong wind conditions. 
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vs. 32, 33, 3^) or wind speed (figures 35 vs. 36) as would be 
qualitatively predicted by a Gaussian model. This model requires 
one additional user input over those needed in the previous 
version; the standard deviation of the wind direction used in the 
meander calculations. As stated in the previous sections, the 
averaging time must be the cloud travel time to the distance of 
interest. Also, the value used should be the average of several 
such standard deviations obtained over a significant fraction of 
the day. The sampling time is not critical , 1.0 - 0.1 Hz would 
be sufficient. 

Code modifications are listed in figures 3 7 - ^ 0 . Figure 37 
shows both the turbulence intensity input and the revised sigma 
parameter table (relative diffusion). Figure 38 shows where the 
meander envelope is calculated while figures 39-^0 list the 
plotting routine. 
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420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

3790 

3800 

3810 

3820 

3830 

3840 

3850 

3860 

3870 

3880 

3890 

3900 

3910 

3920 

3930 

3940 

3950 

3960 



a) 



M e an_w i n d_s p_k t =5 ! K T S 

Mean_u i nd_d i r =275 ! DEG TRUE 

i 

M i x_ 1 ay e r_h t=0 ‘METERS 

! 

! $$$$$ NEW 1985 $$$$$$$$$ 

Tur b_i nt ens= 1 2*P I / 1 80 ! turbulence intensity (radians) input 

! $$$$$$$$$$$$$$$$$$$$$$$1 

| 

DIM Ml_ht_est_meth$C30] 

Ml ht est meth$="** DEFAULT 



b) 



! LOAD ARRAY CONTAINING VALUES OF Ap->Dp 

i 

FOR P0=1 TO 6 ! WHERE P0 = NUMERICAL REFERENCE TO STABILITY CLASSES 

' < 1 = " A " -> 6 = "F " ) 

READ Mat r l x_ap < P0 ) , Mat r i x_bp ( P0 ) , Mat r i x_c p ( P0 ) , Mat r i x_dp ( P0 > 

NEXT P0 



$$*$$$$ NEW 1985 $$$$$$$$$ 



DATA 


0, 0, 0, 0 


i 


-- CAT “ 


DATA 


0.0410, 1 


.00,0.32,0.75 


! -- CAT 


DATA 


0.0410,1 


. 00,0.32,0.70 


! -- CAT 


DATA 


0.021 1 , 1 


. 00, 0. 16, 0. 65 


! -- CAT 


DATA 


0.0211,1 


.00,0. 10,0.62 


! -- CAT 


DATA 


0, 0, 0, 0 


! 


-- CAT “ 



s$stsssssst$$ssss$$$$sts$$ 



(NOT USED OVER THE OCEAN) 

B" 

C" sigmas are for relative diffusion 
D" in lateral direction only 
E" 

> " ( NOT DEFINED OR USED BY THIS MODEL) 



Figure 37. Abbreviated listing of CWHFP showing (a) turbulence 
intensity input and (b) new array of sigma parameters. 
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5710 

5720 

5730 

5740 

5750 

5760 

5770 

5780 

5790 

5800 

5810 

5820 

5830 

5840 

5350 

5860 

5870 

5880 

5890 

5900 

5910 

5920 

5930 

5940 

5950 

5960 

5970 

5980 

5990 

6800 

6010 

6020 

6030 

6040 

6050 

6060 

6070 

6080 

6090 

6100 

6110 

6120 

6130 

6140 

6150 

6160 

6170 

6180 

6190 

6200 



! *********************************************************************** 

Pt_vapor_nom i x : ! CALC MODEL FOR POINT SOURCES OF GASEOUS AGENTS, 

! AND NO REFLECTIVE MIXING LAYER CAP 



! $$$$$ NEW 1985 $$$* 

DIM Meand_x < 1 00 > , Meand_y < 1 00) 

DIM Dat a_set_x < 3 , 1 00 ) , Dat. a_set_y < 3 , 1 00 > ! 100 points vs. 20 pts 

Po i nt s_per_set = 1 00 ! in old version 

i 

M_to_yds=39. 37/36 ! CONVERSION FACTOR USED IN EQUATIONS BELOW, DEFINED 

! AGAIN HERE TO PERMIT TESTING OF PARTIAL PROGRAM 

i 

FOR K=1 TO Num_data_sets 

Par t i al _1 < K ) =$ourc e_s i ze <K ) / CP I *60*Mean_w i nd_sp_m*Dose_val C K ) ) 
Max_rangs_x(K) = (Part i al _1 00 / < Ap < K ) *Cp < K ) > ) * <1 / C Bp < K > + Dp < K > ) ) 

D at a_s e t _y < K , 0 ) = 0 

X_increment =Max_range_x OO/Point s_per_set 

IF X_i nc rement =0 THEN X_i nc re ment = 1 ! ADDED TO PREVENT BLOWUP ON BAD DATA 

i 

FOR J=1 TO Points_per_set 
X = J * X_ increment 
Dat. a_set_x < K , J ) = X 

S i g m a_y x = Ap<K)*X^Bp< K> ! ASSUMES NO INITIAL CLOUD SIZE FOR CONSERVATISM 

S i gma_ zx = Cp C K)*X^DpOO 

Partia1_2 00=2*LOG<Partial_100''<Si gma_yx*S i gma_zx> ) 

IF Part l a1_2(K><0 THEN Par t i al _2 < K > = 0 ! PREVENT POSSIBLE ERROR IN LATER SQR 

CALCULATIONS 

D at a_s e t _y < K , J ) = S i g m a yx*SQR(Partial_2<K>) 

NEXT J 
NEXT K 

i 

! *$$** NEW 1985 $***$$$** 

RAD 

FOR J=1 TO Points_per_set 
X = J*X_ increment 

IF X >Max_range_x < Num_dat a_set s ) *COS ( 2*Turb_i nt ens > THEN GOTO Close 
Meand_y < J)=X*TAN<2*Turb_i nt ens) ! meander envelope is based 

Meand_x(J)=X ! only on turbulence intensity and 

GOTO J100 ! selected to be the 2 sigma value <95*/. ) 

i 

Close: Meand_y C J ) =SQR < Max_r ange_x < Num_dat a_sets) ^2-X^2> ! stop envelope 

Meand_xCJ)=X ! at maximum range and close with partial circle 
J100: NEXT J 

MAT Meand_x=Meand_x*(M_to_yds) 

MAT Meand_y=Meand_y * C M_t o_y ds ) 

i 

! nmmmmiimmi 

Figure 38. Abbreviated listing of CWHFP showing calculation 
of the meander envelope. Points per data set may be adjusted 
for speed or better resolution. 
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8190 ! St$*$ NEW 1985 *$SSS$SS$*$ 

8200 DIM Scaled_x(100),Scaled_y<100> ! LIMITS BASED ON "POINT S_P E R_S E T M 
8210 DIM Cor.tour_l abl_x<3> ,Contour_l ab1_y(3> ! LIMIT BASED ON NUM_D AT A_SE TS 

8220 ! 

8230 FOP Poi nt_num = 0 TO Po i nt s_per_set 1 PO I NT$_PER_SET preset, above 

8240 Sc al e d_x < Po i nt _num ) =Meand_x < Po i nt _num ) *Sc al e_f ac t or ! scale meander 
8250 Sc al e d_y Po i nt _num > =Meand_y < Po i nt_nurn > *Sc al e_f ac t or ! envelope 
8280 NEXT Point_num 
8270 GOSUB Plot_meander 
8280 » 

8290 ! 

8300 1 

8310 FOR Set _num = 1 TO Hun._dat a_se t s ! NUM_DATA_SETS PRESET TO "3", ABOVE 

8320 ! 

8330 FOP Point_num=0 TO Po l nt s_per_s et 1 PO I HTS_PER_SET PRESET ABOVE 

8 340 Seal e d_x (Poi nt_num ) = Dat a_set _x <. Set_num , Po i nt_num ) *Sc a 1 e_f ac tor- 
8350 Scale d _y (Poi nt _num > = Dat a_s e t _y ( Set _n urn , Poi nt _num > *Sc a 1 e_f act or 
8 3 6 0 NEXT Point _n u m 
8 370 GOSUB F'lot_contour 

8380 Cent our_l ab l _x (Se t_num ) = Sc al ed_x C Poi nt s_per_s et ) 

8390 ! 

8400 i ***$**$$* NEW 1985 S$tSS$SSSS$t 

8410 RAD 

8420 R = Dat a_se t_x ( Se t_num , Po i nt s_per_set ) / Dat a_se t_x ( Num_dat a_set s , Po i nts_per 
et > ! dimensionless distance from source [0,13 

8430 0 m e g a= 4 0 * P I / 1 8 0 / T u r b_ i n t e n s ! ARBITRARY frequency of ripple oscillation 

8440 R=Sc al ed_x < Po i nt s_pe r_se t > *TAN < Turb_i nt ens > ! ARBITRARY ripple amplitude 

8450 Ripple=A*SIN(Omega*PI*R> 1 offset for y position 

8480 C ont our_l ab 1 _y ( Se t__num ) = Sc al ed__y < Po i nt s_per _set > +R i pp 1 e ! add offset to 

b e 1 position 

8470 DEG 

s*4 80 1 

8490 ! 

8500 NEXT Set_num 

8510 GOTO Ex i t_dat a_p 1 ot 

8520 1 

8530 1 $$*$$*$ NEW 1985 t SSSSSSSSS SSStSSS ****$?*$$*$$* (NEW SUBROUTINE) 

8540 ! 

8550 P 1 ot _me ander : ! SUB to plot meander envelope 

8560 LINE TYPE 1 

8570 Max_overst r i k e= 1 ! same as "old" version 

8580 X_offset=0 

8590 ! TOP HALF 

86O0 FOR St ri kecount = 0 TO Max_overs t r i k e 
8610 Y_of f se t = . 33*St r l k ec ount 

8620 X_of f set =. 33*St r i k ec ount ! strikecounts thicken line 
8630 MOVE Sc al ed_x ( 0 ) , Sc al ed_y ( 0 ) 

8640 FOR K = 0 TO Poi nts_per_set 

8650 DRAW Sc a 1 ed_x ( K > + X_of f set , Sc al ed_y ( K > ♦ Y_of f se t 
8660 NEXT K 

8670 ! SCALED_X,Y now contain meander envelope 

8680 i PLOT "BOTTOM" HALF 

8690 MOVE Seal ed_x<0> , Seal ed_y(0> 

8700 FOR K = 0 TO Poi nts_per__sct 

8710 DRAW Seal ed_x(K>+X_of fset , - ( Sc a 1 ed_y OO ♦ Y_of f set ) 

8720 NEXT K 

8730 NEXT Strikecount 
8740 RETURN 

8750 ! $$$$$$$$$$$$$$$$$$$$$$$$$$$ 

Figure 39. Abbreviated listing of CWHFP showing meander envelope 
plotting routine and contour label locating scheme. 
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8770 ' ******* NEW 1985 ********** (REVISED SUBROUTINE; 

8780 PI ot_contour: ! plot, contours of the scaled data 

8790 RAD 

8800 LINE TYPE 1 

881© Gniega=4u*P I / 1 80 Turb_i nt ens ! ARBITRARY frequency of ripple oscillation 
8820 ! 

8830 ! PLOT BEGINNING WITH "TOP" HALF OF CONTOUR 
8840 ! ( OVERSTRIKE STUFF OMITTED > 

8850 MOVE Sc al ed_x < 0 > , Sc al ed_y ( 0 > 

8860 FOR K=0 TO Poi nts_per_set 

8870 R = Dat ( Set_num , K > •' Bat a_se t_x ( Num_dat a_set s , Po i nt s_per_set ) ! d i mens i 

on less distance from source 

S8S0 A = $c a 1 ed_x OO *T AN < Tur b_i nt en« > ! ARBITRARY amplitude of ripple selected to 
be 1 sigma 

8890 Ripple=A*SIN<Omega*PI*R> ! offset in y direction 
8900 DRAW Seal ed_x<K> , Seal ed_yOO+Ri ppl e ! add offset 
8910 NEXT K 
8920 ! 

8930 ! PLOT "BOTTOM" HALF most same as "top" half 

8940 ! 

8950 MOVE Seal ed_x (0) , Seal ed_y(0) 

8960 FOR K=0 TO Po i nt s_per _se t 

8970 R = Dat. a_set _x ( Set _n urn , K > / Dat a_se t_x ( Num_dat a_set s , Po i nt s_per set ; 

8980 A = Sc al ed_x <K>*TAN<Turb_intens) 

8990 Ri pple=A*SIN(Omega*PI*R) 

9000 DRAW Sc al ed_x < K ) , - < Sc al ed_y < K > -R i pp 1 e > ! subtract offset 

901O NEXT K 
9020 DEG 

903O RETURN ! END OF GOSUB PL0T_C0NT0UR 
9040 * ***$**$$***$*********$***$*$**$***** 



Figure 40. Abbreviated listing of CWHFP showing plotting of 
hazard contours based on relative diffusion parameterization. 
"Wiggles” are a result of a SIN wave imposed on the Gaussian 
plume model solution. 
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6. SUMMARY AND FUTURE WORK 



The presented work successfully parameterizes relative 
lateral diffusion from mean meteorological quantities, verified 
via two independent overwater tracer data sets. Single-particle 
diffusion is successfully parameterized only through direct 
measurements of the lateral turbulence intensity. Once this 
quantity is known, however, the parameterization is well behaved. 

Future work should attempt to correlate this lateral 
turbulence with large scale synoptic features, mesoscale 
phenomena, radiosondes profiles, or other more easily measured 
quantities. The GMGO researchers are presently being petitioned 
for positioning information on their experiments so that a more 
thorough verification of the single-particle parameterization can 
proceed in the future. 

The model is presently designed for only medium ranges, but 
could be extended to greater distances. This extension would 
make results highly dependent on inversion height. Some 
prognostic or diagnostic estimate of this quantity could be 
incorporated into the model. Longer range overwater experiments 
should be researched for plume paramterizations at these 
distances . 
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